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Abstract 

We study models of gauge mediation with strongly coupled hidden sectors, employing 
a hard wall background as an holographic dual description. The structure of the soft 
spectrum depends crucially on the boundary conditions one imposes on bulk fields at the 
IR wall. Generically, vector and fermion correlators have poles at zero momentum, leading 
to gauge mediation by massive vector messengers and/or generating Dirac gaugino masses. 
Instead, non-generic choices of boundary conditions let one cover all of GGM parameter 
space. Enriching the background with R-symmetry breaking scalars, the SSM soft term 
structure becomes more constrained and similar to previously studied top-down models, 
while retaining the more analytic control the present bottom-up approach offers. 



1 Introduction 



We have started in [T] to investigate with holographic methods [2j|7] strongly coupled 
hidden sectors with the aim of extracting data suitable for models of gauge mediation of 
supersymmetry breaking, using the formalism of General Gauge Mediation (GGM) [8]. 
We have laid out the general procedure one should follow, and applied it to a class of 
non-supersymmetric Asymptotically AdS (AAdS) backgrounds which arise as solutions 
of ten-dimensional type IIB supergravity (9j[l0). As such, these backgrounds represent 
a top-down approach where the supersymmetry breaking theory is found directly in a 



possibly string-derived and UV-complete set up (see 11-14 for other works in the same 
direction) . 

The two rather general results we have obtained there can be summarized as follows: 

• Pure dilatonic backgrounds describe, holographically, supersymmetry breaking hid- 
den sectors where gauginos cannot acquire a Majorana mass. Rather, they can get 
Dirac masses, due to direct coupling to R-charged massless fermions, which arise 
quite generically in the low energy spectrum of the strongly coupled hidden sector. 
As such, the supersymmetric Standard Model (SSM) soft spectrum turns out to be 
similar to that of gaugino mediation models. 

• Backgrounds where scalars whose field theory dual operators are charged under the 
R-symmetry get a non-trivial profile, do instead generate gaugino Majorana masses. 
In such situations, no R-charged massless fermions show-up in the hidden sector, 
and hence there is no Dirac-like contribution to gaugino masses. Here the SSM soft 
spectrum looks democratic between gauginos and sfermions, and hence similar to 
that of minimal gauge mediation models. 

A natural question one might ask is whether these results are generic for models 
admitting holographic duals, or whether they strongly depend on the specific background 
we have been considering. More generally, how large a portion of the GGM parameter 
space can holographic models of gauge mediation cover? Are there any restrictions? 
These questions are difficult to answer within fully-fledged top-dow models, given also the 
poor number of concrete and sufficiently explicit string theory supersymmetry breaking 
solutions available in the literature. 

In the present paper we will focus on a complementary, bottom-up approach. We 
will consider simple backgrounds which do not have necessarily a completion in string 
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theory, but which have the advantage to let one have more flexibility and more analytical 
power. This will enable us to answer the above questions, postponing to future studies the 
possibility that the extra requirement of a full string embedding might eventually restrict 
further the region of GGM parameter space one can cover via holographic models. 

We will keep considering AAdS backgrounds, and, without loss of generality, focus 
on a hidden sector global U(l) symmetry. This is the symmetry one should weakly 
gauge in order to couple the (supersymmetry breaking) hidden sector dynamics to the 
SSM. We will focus on a simple hard wall model 15 17] which consists of a pure AdS 



background with an IR cut-off in the bulk. The boundary conditions one imposes on the 
fluctuating fields at the IR cut-off are the only freedom one has in this set up. We will 
parameterize the structure of the SSM soft spectrum accordingly. A scan of all possible 
boundary conditions, which can be parameterized in terms of a number of (a priori free) 
parameters, gives the following results: 

• Generically, that is if one does not tune any of the parameters, the resulting soft 
spectrum is that of gauge messengers mediation scenario. Gauginos acquire Dirac 
masses, while sfermions are loop suppressed, hence the soft spectrum turns out to 
be similar to that of gaugino mediation models. These features arise due to the 
presence of massless (composite) bosonic and fermionic modes in the hidden sector. 

There are two possible ways out of this scenario, whose effects we study in detail: 

• Tuning some of the parameters one can eliminate hidden sector massless modes: 
in this case one can scan any value of sfermion and gaugino masses, and obtain a 
spectrum which can range from gaugino mediation to minimal gauge mediation or 
scenarios with suppressed gaugino masses, hence covering all of GGM parameter 
space. 

• A more constrained but attractive possibility, relies on adding dynamical degrees 
of freedom in the bulk. Within this approach, the hidden sector massless fermionic 
modes disappear automatically, and standard Majorana masses for gauginos are 
generated, without any necessary tuning of bulk field boundary conditions. 

The basic objects one has to compute in GGM are two-point functions of operators 
belonging to a hidden sector conserved current supermultiplet. For a £7(1) current these 
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are parametrized in Euclidean signature and in momentum space as 



(J(k)J(-k)) = C (k 2 ) , 
(Uk)M-k)} = -vLkiC 1/2 (k 2 ) , 
(ji{k)jj{-k)) = {kikj - Si^C^k 2 ) 
UMjp(-k)) = e a pB 1/2 (k 2 ) , 



;i.2) 

1.3) 

;i-4) 



where a factor of (27r) 4 5^ 4 ^(0) is understood. Gaugino and sfermion masses can be ex- 
pressed in terms of the above correlation functions as[j] 



m x = g 2 B 1/2 (0) 



171: 



2 2 

9 Qf 



E 



\2.s 



N„ 



1.5) 
1.6) 



where qj is the sfermion U(l) charge, N = 1, Ni/ 2 — 4, Ni — 3, and g is the U(l) gauge 
coupling. We have written the expressions ( 1.5 )— ( 1.6 ) following the spirit of [l8j[l9], in 
order to allow for massless poles in the C s correlators. We will indeed discuss cases where 
some or even all of them will have such a singularity. When all correlators are instead 



non-singular in the IR, one can expand the expression (1.6 ) perturbatively in g, and obtain 



the usual leading order GGM expression for the sfermion masses derived in [8]. 

The program of holographic general gauge mediation outlined in [lj consists in com- 



puting the correlators ( 1.1 )-( 1.4 ) in strongly coupled hidden sectors using holography. 



1.1 Outline 

We begin in section 2 by recalling the techniques of holographic renormalization for AAdS 
spaces, when applied to a massless vector bulk supermultiplet. The latter is the holo- 
graphic dual of the conserved current supermultiplets whose two-point functions we want 
to compute. This will allow us to provide an holographic expression for the two-point 



functions (1.1)-(1.4), for a generic AAdS background. In order to set up notations, and 
as a useful reference for later studies, we will then apply such general formulas to the 
(superconformal) pure AdS case. In section 3 we move to consider the hard wall model, 
and show that one can impose a variety of different boundary conditions on the bulk 



Strictly speaking, eq.(1.5) is approximate, the exact expression being given by the pole of the re- 
summed gaugino propagator. The exact value can differ significantly from the approximate one when 
g 2 Bi/2 and/or g 2 Ci/2 are 0(1), as in some of the strongly coupled (large N) scenarios that we consider. 
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fields at the IR wall, obtaining a corresponding hidden sector in which the breaking of 
both supersymmetry and R-symmetry depends on independent parameters in the bound- 
ary conditions. The phenomenological consequences of this freedom will be discussed in 
detail in section 4, where it will be shown when and how the general results we have 
anticipated above, can actually be achieved. A somewhat unpleasant feature of this way 
of proceeding is that in order to cover large portions of GGM parameter space, and in 
fact allow for non-vanishing Majorana gaugino masses to begin with, somewhat ad hoc 
boundary conditions should be imposed. In section 5, by turning on a background profile 
for a scalar whose dual operator is charged under the R-symmetry, we will show that 
gaugino masses can in fact be generated in a dynamical way, following a similar strategy 
as that of [l] . We conclude in section 6 with a brief summary and an outlook on further 
research directions. The appendices contain a list of useful properties of Bessel functions 
that we use in our analysis, and the derivation of some analytic results about the behavior 
of GGM correlators at low momenta which we will use in the main text. 



2 Holographic correlators for GGM 

Let us summarize the main requirements that a background has to satisfy in order to 
be interesting from our perspective. These requirements are directly related to specific 
features of the hidden sector gauge theory that we want to describe holographically: 

1. The requirement of Asymptotically Anti-de Sitter-ness results in having a gauge 
theory with a superconformal (interacting) UV fixed point. This is not a necessary 
physical requirement, but just a simplifying assumption we make to have more 
control on the AdS/CFT machinery, and which can be relaxed if one wants to 
consider other classes of (more realistic) hidden sectors. 

2. The background should asymptote AdS but deviate from it in the interior. This 
ensures the existence of one or more scales in the problem, which is a necessary con- 
dition for the background to break conformal invariance and supersymmetry. Single 
scale models are the simplest option, but more realistic gravitational duals might 
admit different (and unrelated) scales setting the breaking of conformal invariance 
and supersymmetry, respectively. Typically, the field theory will enjoy a confining 
behavior, with a mass gap. 
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3. The background has to support an iV = 2 vector multiplet corresponding to a 
gauged symmetry group - at least an U(l) - in the bulk. This corresponds to an 
anomaly-free global symmetry of the hidden sector, a necessary ingredient for a 
model of gauge mediation, where this global symmetry should eventually be weakly 
gauged. 



The starting point for computing the GGM correlators ( 1.1 )-( 1.4) holographically is 
to turn on the fluctuations for the fields dual to the conserved current supermultiplet, 
which, as just noticed, form an N = 2 vector multiplet in AdSs. Since we are interested 
in two-point functions, we need the action for the fluctuating fields only to quadratic 
order. In what follows, we will start considering a minimal action, with kinetic plus mass 
terms only 

N 2 



Smin = J I d XVG 



{GTdpDd v D - AD 2 ) + l -F^ + 1 -{\H) G \ + c.c.) - ^AA 



(2.1) 

where we put the AdS radius to be L — 1, and the masses = 0, m 2 D = —A and 
m\ = —1/2 are those required by the holographic dictionary. The equations of motion 
are 

(□ G + 4)D = 0, (2.2) 

-^=d^VGGWG™F pa ) = , (2.3) 
V G 

{U) G - = . (2.4) 

In a general AAdS background the asymptotics of the fluctuating fields in the vector 
multiplet is uniquely determined by their masses. In the coordinates 

dsl = \ (dz 2 + dx 2 ) , (2.5) 

the boundary is at z — > 0, and near the boundary the equations of motion become 

(z 2 d 2 - 3zd z - z 2 k 2 + A)D = , (2.6) 
(z 2 d 2 - zd z - z 2 k 2 )A t = , (2.7) 



[z 2 d 2 - Azd z - z 2 k 2 + f )£ = 



za l kiX = (zd z - §)£ 



4 

3^ 



(2.8) 
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In order to obtain the equations above we have Fourier transformed the fields, and we 
have fixed the gauge for the 5d gauge field taking A z = and —ik l Ai = 0. Moreover, 
we have traded the first order equation of motion for a Dirac field with a second order 
equation of motion for one of its Weyl components plus a first order constraint for the 
other Weyl component, with 

The resulting near boundary expansion is 



(2.9) 



D(z, k) ~ z 2 (\n(eA)d (k) + do(fc)) + 0{z 4 

z— s-0 



A- 2 



Ai(z,k) ~ a i0 (k) + z 2 (a i2 (k) + — \n(zA)a i0 (k)) +0(z 4 ) , 
z->o 2 

~ z V\^k) + z\i 2 (k) + k 4HzAMk)) + 0(z^) 
X (z,k) ~ ^/ 2 (xi-ln(zA)a J ^o(A;) + 0(^)) . 



(2.10) 
(2.11) 

(2.12) 



The leading terms in the expansion {do(k), aio(k), £q (k)} will be identified as the sources 
for the corresponding boundary operators {J(k), ji(k), ja(k)} and the undetermined sub- 
leading terms {do(k), ^(k), Xia{k)} will be related to the one-point functions of the same 



boundary operators. The near boundary limit of the action (2.1 ) appropriately renormal- 



ized gives us the 4d action from which one can derive the holographic expression for the 
GGM two point functions 



(J(k)J(-k)) = 

(ja(k)ja(-k)} 
(ja(k)M-k)) 




2 ff + k% 

6a l 



(2.13) 
(2.14) 



Id 



8^ 



1/3 




k 2 



(a <H- (3) 



(2.15) 
(2.16) 



The fermionic correlators (2.15) and (2.16) are naturally expressed in terms of the sublead- 



ing mode Xi(k) of the positive chirality field x{ z i k). For later convenience, we have also 
re-expressed the formulas in terms of the subleading mode &(k) of the negative chirality 
field, which is related to Xi(fc) by the equations of motion. 
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2.1 Pure AdS: the superconformal reference results 



As a reference for what we do next, let us now review the standard supersymmetric AdS 
computation of [j] . The equations of motion for the fluctuating fields in AdSs are exactly 



the second order linear ODE (2.6)-(2.8) defined in the open set (0,oo). Performing a 
rescaling of the bulk fields we can bring these equations to the standard form for the 
Bessel equation (see appendix [A] for a review of the essential features of Bessel functions). 
In fact, taking 

zai , I = z^E , X = z^X , 



D = z 2 d , Ai = 
we get for the equations of motion 

{z 2 d 2 z + zd, 

(z 2 dl + zd. 



z l k l )d = , 
(z 2 k 2 + i))oi = 



{z 2 d 2 + zd z - {z 2 k 2 + 1))S = , 
za%X = (zd z + 1)2 , 

whose general solutions for k 2 > are 

d{z,k) = ci(k)I (kz) + c 2 (k)K (kz) , 
aci(z, k) = a a (k)Ii(kz) + a i2 (k)Ki(kz) 
E(z, k) = B x {k)h{kz) + B 2 {k)K x {kz) . 



(2.17) 

(2.18) 
(2.19) 

(2.20) 



(2.21) 
(2.22) 
(2.23) 



The general solutions of second order differential equations depend on two arbitrary con- 
stants, that in our case can be arbitrary functions of the Ad momentum k. 

In order to impose Dirichlet boundary conditions at z — » 0, which identify the sources 



for the boundary operators, we expand the solutions (2.21 )-(2.23) using (A.9)-(A.10) 



Comparing with (2.10)-(2.12) we get 



c 2 (k) = -do (&) , oi i2 = ka i0 (k) , 9 2 = k£ (k) . 



(2.24) 



This leaves us with three arbitrary functions c\(k) , an(k) , 6i(k). In order to have a regular 
solution in the full domain we have to impose the following conditions in the deep interior 



lim D(z, k) = , lim A^z, k) = , lim £(z, k) = . 



(2.25) 



Using the expansions (A.7)-(A.8) it is easy to see that (2.25) implies c\{k) = an(k) 



6i(k) = 0. We note that the bulk boundary conditions are crucial to single out the solution 
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to the fluctuation equations. The pure AdSs solutions are 



D AdS (z,k) = -z 2 K (kz)d (k) 
A AdS (z,k) = zkK x {kz)a iQ {k) , 
Z AdS (z, k) 
AdS {z } k) 



X 



z"' 2 kK x {kz)Uk) , 
z 5 / 2 K (kz)a%Co(k) 



(2.26) 
(2.27) 
(2.28) 
(2.29) 



where in order to derive the last equation we have used the recurrence relations among 



the Bessel functions (A.6). 



Expanding these results near the boundary we get the usual AdS 5 result for the sub- 
leading modes from which we derive the GGM two-point functions 



(J(k)J(-k)) = 

(mjA-k)) - 

(ja(k)]a(-k)) 

(Uk)jp(-k)) 



N 2 



a 



AdS, 



8tt 2 ^ 
N 2 



in 
N 2 



(k 

2$ijk Ci ^ (k 



AT 2 

2 



87T 

2 r< AdS ii 2 



In 



57T 



2 a%C A f(k 2 ) 



A* 

N 2 
8^ 
N 2 



+ 21n2-27 



2 Sijk 2 



cr'ki 



In 



In 



A? 

k? 
A 2 
k 2 



+ 2 In 2 - 27 
+ 21n2 -27 



(2.30) 

(2.31) 

(2.32) 
(2.33) 



where we can always take in (2.31) 5ij — > 5ij — j^- because of current conservation^ 
As expected in a supersymmetric background, Bi/ 2 = and CQ dS = C AdS = C Ad 2 s 



C AdS so that both the Majorana gaugino mass (1.5) and the sfermion masses (1.6) are 



identically zero. It is worth noticing that the correlators (2.30)-(2.33) reflect the fact 
that the hidden sector is exactly superconformal, exhibiting a branch-cut at k = which 
corresponds to the two-particle exchange of the 2A^ 2 massless particles of the hidden sector 
which are charged under the gauged U(l). 



3 Hard wall models 

Let us now consider the model that will be our primary interest in this paper, the hard 
wall model. This is just AdSs in which the geometry ends abruptly in the interior by 
putting a sharp IR cut-off at z = This model was originally studied as a toy model 
of a confining gauge theory because it provides an holographic dual for theories with a 

2 As was already noticed in [l] the holographic prescription performed with our gauge fixing allows us 
to compute only the part of the vector current two-point function which is proportional to 8ij. 
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gapped and discrete spectrum |15f|17| . Our aim is to study the behavior of GGM two- 
point functions on this background]^] Not surprisingly, the behavior of the correlators will 
depend strongly on the boundary conditions that one has to impose on the bulk fields at 
the IR cut-off. 

In the case of an hard wall background the general solution of the equations of motion 
for the fluctuations in the vector multiplet is exactly the same as for pure AdS, eqs. 



(2.21)-(2.23), and depends on six integration constants (two for each field). Also the UV 



boundary conditions (2.24) remain the same, and can be simply understood as fixing the 
source of the boundary operator, leaving only three constants undetermined. However, 
we are now solving the differential equations in the open set (0, 1/V), and the regularity 
conditions are replaced by some IR boundary conditions at z — 1///. These conditions 
can be solved for the three remaining constants, in order to fix the functional dependence 
of the subleading modes on the leading ones. 
The generic form of the solutions is 

D(z, k) = D AdS {z, k) + Cl (k)z 2 I {kz) , (3.1) 
A(z, k) = A AdS (z, k) + a u {k)zh{kz) , (3.2) 
k) = l AdS {z, k) + e^z^hikz) , (3.3) 

where ci, an and d\ are integration constants determined by the IR boundary condi- 
tions. Expanding these expressions near the UV boundary, one can easily find that the 
correlators are modified with respect to the AdS case in the following way 

C (e)=C AdS (k*)-2 d £ , (3.4) 

0l(O =C-V)-l*,g, (3.5) 

C 1/2 (fc 2 ) = C" s (k>) - f^f |f + C A • < 3 - 6 ' 

B^ 2 ) = -^ir- (3-7) 
k b£ 0a 



Let us notice that our set up is reminiscent of extra dimensional scenarios like 20 



22 



in which, 



however, the physics of the Ad hidden sector arises as a KK reduction of a 5d theory in a slice of AdS. 
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3.1 Homogeneous IR boundary conditions 

We start by taking general homogeneous boundary conditions at the IR cut-off 



(D(z, k) + Po zd z D(z, k))\ z= i,p = , 
{Ai(z, k) + p 1 zd z A i (z, k))\ z=1/lx = 
k) + pi/ 2 zd z £(z, k))\ z=1/fl = , 



(3.8) 
(3.9) 
(3.10) 



which depend on three independent coefficients p s . As we will see, in order to cover all of 
GGM parameter space it will be necessary to turn on also inhomogeneous terms in the 
above equations, something we will do next. 

As it befits coefficients computed with homogeneous boundary conditions, the coeffi- 



cients Ci, an and d\ in (3.1)-(3.3) are all proportional to the source terms. The resulting 
GGM functions are 



c[ h \k 2 ) 

(k 2 ) 



+ 2 



'1/2 V 



a 

B[%(k 



C AdS (k 2 ) 
C AdS (k 2 ) 
C AdS (k 2 ) + 2 



-(l + 2p )K {±)+p ±K 1 ( k 



2po)/o(J)+A)*/i(5: 



+ 2- V 



(l + |pi/2)ATi(g)-Pi/2^o 

(l + |Pl/2Kl(^)+Pl/2^0( 



. 



(3.11) 

(3.12) 

(3.13) 
(3.14) 



The analysis of the boundary condition-dependent soft spectrum emerging from the 
correlators (3.11 )— (3.13) is postponed to section |4j For future reference we would instead 
like to comment here on the behavior of the correlators in the IR and UV. Making use 



of the asymptotic expansion for i C 1 for the Bessel functions (A.9)-(A.10) we find the 
correlators at low momentum to behave as 

'A 2 ' 

72 



C ( h \k 2 



~ In 

k^O \ p' 



+ 



2p 



C[ h \k 2 



C[%(k 2 ) 



k^O 1 



2pi k 



l + 2p 



2 + 3p 1/2 p 2 
~ 4 '—!—'— -fin 

fc->o 2 + 7pi/2 k 2 



_ 3 + 8p! 
"2(l + 2 Pl )2 ' 
A_\ (2 + 3p 1/2 )(6 + 25p 1/2 ) 
A*V 2(2 + 7 Pl/2 )2 



(3.15) 
(3.16) 
(3.17) 



As for the UV limit, given the large x behavior of Bessel functions (A.7)-(A.8), we can 



see that all the functions approach the supersymmetric AdS value with exponential 
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rate at large momentum 

C^ h \k 2 ) ~ C[f 2 (k 2 ) ~ C[ h \k 2 ) ~ C Ad5 (A; 2 )-27r e _ ^. (3.18) 

From the field theory point of view, the exponential suppression of the SUSY-breaking 
effects in the UV suggests that SUSY-breaking in a hidden sector described by a hard 
wall holographic model is not induced by any quantum operator, which generically would 
appear in the OPE of the GGM correlation functions with a scaling behavior in k 2 fixed 
by its dimension. 

Two additional remarks are in order at this point. The first is that one can of course 
compute the above functions also using the numerical approach pursued in fij, finding 
perfect agreement with the analytic computation above. The second comment is that 
the above functions can be continued to negative values of k 2 . It is easy to convince 
oneself that they will then display an infinite sequence of poles on the negative k 2 axis, 
corresponding to the glueball towers for each spin sector. They return the same values 
that can be obtained through the more traditional holographic approach of computing 
glueball masses, i.e. finding normalizable fluctuations for each field. 



3.2 Inhomogeneous IR boundary conditions 

Let us now consider the possibility of having inhomogeneous boundary conditions in the 
IR. We thus take general boundary conditions at the IR cut-off depending on three more 
arbitrary terms, now 

(D(z, k) + Po zd z D(z, k))\ z=1/fl = E (k) , (3.19) 
(Ai(z, k) + pizd z Ai{z, k))\ z=1/fl = E n (k) , (3.20) 
(£(z,k) + pi/2zd g £(z,k))\ z=1/lt = £ 1/2 (£;) , (3.21) 

where we have allowed for a non trivial k dependence in the inhomogeneous terms E s . 
We will see instantly that the arbitrariness actually amounts to 4 new constants. 



The coefficients c\, an and 9\ in eqs. (3.1)-(3.3) will pick up an additional contri- 



bution, linear in the E s . Since these coefficients enter the GGM functions only through 
the first derivative with respect to the source, the inhomogeneous terms can contribute 
only if we allow them to be dependent on the source, with the result that the condition 



at z — 1/p: involves both IR and UV data of the function. In particular, from eq. (3.7), a 
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dependence of S 1 / 2 (A;) on the source £o can give a non-vanishing By 2 , as opposed to the 



case of homogenous boundary conditions (3.14). Therefore, such a dependence implies 



that the boundary condition (3.21) explicitly breaks the R-symmetry. 



Since in any case only the first derivative enters eqs.(3.4)-(3.7), it is enough to let the 
S s depend linearly on the sources do(k) , aio(k) and £o(^)- Taking into account Lorentz 
invariance, a reasonable choice is 

S (A; 2 ) = --EMk) , (3.22) 
p 2 

^(k 2 ) = -E 1 a iQ {k) , (3.23) 

^iMk 2 ) = -~E 1/2 $(k) - Hy^ark'Uik) , (3.24) 

where the £"s and H are coefficients which do not depend on the momentum. Hence we 
are left with 4 new parameters due to the inhomogeneous boundary conditions. 
The GGM functions in this case take the form 

C^k 2 ) = W ) + 2^ + 2 ^:^ + E ° , (3.25) 
C^{k 2 ) = C AdS {k 2 ) + 2 ^%^^||^ , (3.26) 

C (n h)(k2) _ cAdS(k2) (l + lPl/2)^(^)-Pl/2^0(g) + ^l/ 2 

1/2 ( ) ( * ) + 2 d + lPv^ + ^e) ' (3 ' 27) 

< )( " 2) = 2 fl + 3, kl(k) • (3-28) 

The result with homogeneous boundary condition is simply recovered by setting the E's 
and H to zero. 

The inhomogeneous terms contribute to the IR behavior as follows 

Ct h \k 2 ) - C^\k 2 ) c —^—E , (3.29) 
C^-C^^^p^ (3.30) 
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In particular, having Hi/ 2 7^ we get now a non-zero Majorana mass for the gaugino. 
Indeed, the boundary condition (3.24) explicitly breaks the R-symmetry. 



As for the UV asymptotic, the large x behavior of the Bessel functions (A. 7) tells 
us that the exponential approach to the supersymmetric limit remains valid in this case, 
also for Bi/ 2 {k 2 ) that asymptotes to 0. So we see that, consistently, the inhomogeneous 
boundary conditions do not modify the UV behavior. 



4 Analysis of the soft spectrum 

We now discuss the physical interpretation, in terms of soft supersymmetry breaking 
masses, of the C s and B functions we have found in the previous section. 

Let us start with a very basic requirement: since the correlators happen to have non- 
trivial denominators which depend on the momentum, we should exclude the possibility 
that tachyonic poles are developed. The denominators are linear combinations of two 
Bessel functions evaluated at x = k/fi, and studying their monotonicity properties and 
their limits for x — >■ and x — > 00 one can easily see that the poles are excluded if and 
only if the coefficients of the linear combination have the same sign. This condition results 
in the following inequalities 

{Po < -\] U {p > 0} , { Pl > 0} , {p 1/2 < —} U {p 1/2 > 0} . (4.1) 



The IR behavior of the C s functions, in particular the expressions given in eqs. (3.15) 



(3.17), shows that the theory described holographically by the hard wall has a threshold 
\i for the production of two particle states and possibly a certain number of massless poles 
which depends on the choice of the boundary conditions. Below we analyze the cases of 
homogeneous and inhomogeneous boundary conditions in turn. 

4.1 Homogeneous boundary conditions 



For generic choices of p s parameters, we see from eqs. (3.15)— (3.17) that C\ and C1/2 have 
poles at k 2 = while Cq has not. The interpretation of such poles is that they arise from 
the exchange of a massless state with the same quantum numbers of the corresponding 
operator. In C\, this means that the global U(l) symmetry is spontaneously broken, the 
massless excitation being the associated Goldstone boson. If the symmetry is broken, we 
cannot identify it with the Standard Model gauge group, but rather with an extension 
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thereof by some higgsed U(l)', a setting extensively studied in the literature (see for 
instance [23 J and references therein). 

The pole in Cu% signals the existence of an R-charged massless fermion, neutral under 
the global £7(1), which mixes with the fermionic partner of the current. The most natural 
interpretation of such a fermion in a strongly coupled theory is that of a 't Hooft fermion 
associated with a global anomaly of the unbroken U(1)r. 

The consequence on the soft spectrum of poles in the correlators C\ and/or Ci/ 2 
was studied in 18,19 , and can be summarized as follows: the gaugino acquires a Dirac 
mass by mixing with the would-be massless fermion in C\i 2 (recall that a Majorana 
mass is forbidden by the unbroken R-symmetry), and the integral giving the sfermion 
masses is dominated by the contribution of the poles. Comparing with the usual result 
in General Gauge Mediation without IR singularities, the sfermion soft mass is enhanced 
by a logarithm of the gauge coupling. Notice that the pole in C\ (C1/2) contributes with 
a negative (positive) sign, so that generically one can get a tachyonic contribution to the 
sfermion mass-squared. In formulae 



m 2 j: ~ 



gM 1/2 
(4tt) 2 



r 



log^ ) (4Mf /2 -3Mt 



(4.2) 
(4.3) 



where g is the gauge coupling, rrig is the Dirac mass of the gaugino, mj is the sfermion 
mass, and M 2 is the residue of the massless pole C s ~ M 2 /k 2 . From eqs. (3.16)— (3.17) 
we see that in our mode0 



Ml = 4/i 2 



1 



l + 2pi 



M 2 1/2 



4 M 



2 l + fPl/2 
1 + fPl/2' 



(4.4) 



Notice that in the tachyon-free range (4.1) the two residues are always positive. If we 



further impose the contribution to the sfermion mass-squared (4.3) to be positive, we get 
the additional inequality 

*>4^. (4-5) 

2^1/2 



1 



We see from eqs. (4.2)-(4.3) that in this scenario the sfermions are somewhat lighter that 



the gaugino. This is typical of Dirac gaugino scenarios 24 , though in our model the 
Dirac partner of the gaugino is a strongly coupled composite fermion. 

4 Here and in the following we tacitly assume that the prefactor N 2 /8ir 2 can be set to unity. For this 



value the overall normalization of the physical correlation functions (l.l)-( 1.4) coincides with the one 



that we have used throughout ( 2.13 1— ( 2.16 ) . 
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Tuning the p s parameters 



We now briefly mention different possibilities to evade the generic scenario presented 
above, which can be realized by choosing specific values for the p s parameters. 

1. As a first possibility, consider the case in which M\ 

Pl/2 



ML that is 



Pi = 7~~i~ , (4-6) 
1 + 2P1/2 

while po is kept generic. We are still in a scenario in which the global symmetry is 
spontaneously broken in the hidden sector, and the soft spectrum is described by 
the same formulae as before (notice however that the contribution to the sfermion 
mass-squared is positive, now). Nevertheless, in this case we can argue a different 
interpretation of the physics in the hidden sector, the reason stemming from a 



somehow surprising fact: the condition (4.6) that makes the two residues coincide, 
actually renders the whole C\ and C\/2 functions (3.12 ) and (3.13 ) equal for all values 



of k 2 . As a consequence, one is led to interpret the massless fermion as the partner 
of the Goldstone boson associated to the broken global symmetry, rather than a 't 
Hooft fermion. Since Co differs from C\ = C1/2 for generic po, supersymmetry is still 
broken in the hidden sector, but mildly enough so not to lift the fermionic partner 
of the Goldstone boson. 

2. As a subcase of 1, consider in addition to tune the po parameter to po = —1/2. In 



this case the low momentum expansion (3.15) is not valid, and by repeating the 
analysis one finds that also Co develops a 1/k 2 pole, with residue M5 2 = 4p 2 . As 



explained in 18 19 , a pole in Co is unphysical, unless the hidden sector breaks 
the global symmetry in a supersymmetric manner, so that Co = C1/2 = C\ and a 
massless Goldstone mode is present in all three functions^} Indeed, if we require 
Ml = M\ j2 = Ml, that is p x = p 1/2 = 0, we find from eqs. (|3.11|) 13- 13h that this 
condition is sufficient to ensure Co = Ci/ 2 = C\ for all values of k 2 , supporting the 
interpretation of a supersymmetric global symmetry breaking in the hidden sector. 

3. Finally, p\ and pi/ 2 can also be (independently) tuned in such a way to eliminate 
the massless pole in C\ and Ci/ 2 respectively, the specific values being p\ = expand 



5 In the simple example of a U(l) broken by the VEV of a charged chiral superfield the pole in Cq is 
related to the modulus of the complex scalar. 

6 A global parametrization which avoids infinities could be conveniently given in terms of angles a s , 
the change of variable being p s = tg(a s ). 
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Pi/2 = —2/3. If only one of the two parameters is tuned, the soft masses and the 
interpretation of the physics in the hidden sector remains the same as in the previous 
section, with the only difference that M\ or My 2 are tuned to 0. It is therefore 
more interesting to consider the possibility that both parameters are tuned: in 
this case none of the C s has an IR singularity and we are in a situation similar 
to ordinary GGM, as far as sfermion masses are concerned (the gaugino remains 
massless because the hidden sector does not break the R-symmetry). Since at large 
k all the C s approach their supersymmetric value exponentially, the weighted sum 
— (Co — 4Ci/ 2 + 3Ci/ 2 ) goes to zero at the same rate, so that we can determine the 



sign of the sfermion mass-squared by studying its IR limit. From eqs. (3.15)— (3.17) 



we see that the leading term, with the present values of pi and p 1 / 2 , is given by 

-(C0-4CV2 + 3CV) ~ -t^-- (4.7) 

fc->o 1 + 2p 



In the tachyon-free range (4.1) this expression is negative. Therefore, in this tuned 
scenario we find vanishing gaugino mass and tachyonic sfermion mass. We will 
see later that both this unwanted features can be overcome: one way, which is 
somehow more ad-hoc, consists in enlarging the parameter space by considering 
inhomogeneous boundary conditions; the other, which is more dynamical, consists 
in turning on a R-breaking scalar on top of the hard wall background. Most of what 
follows will therefore consist in improvements of this setting with tuned pi and pi/2- 

4.2 Inhomogeneous boundary conditions 



Let us proceed by considering the functions (3.25)-(3.27), which we obtained adding 
source-dependent inhomogeneous terms in the boundary condition. Besides the p s , we 
have now four additional real parameters to play with, namely the dimensionless E s and 
the R-breaking parameter Hi/ 2 , which has dimension of a mass. 

For generic values of the parameters the situation is analogous to the one with homoge- 
neous boundary conditions, so that the E s parameters appear to be somehow redundant: 
C\ and C1/2 have a massless pole, while Co has not. The major difference with respect to 
the previously considered case is that now H\ii gives a non-zero Majorana mass to the 
gaugino, 

m~ 9 = ^-4—^i/2 • (4.8) 



16 



Since now R-symmetry is broken, the pole in C\i% cannot be interpreted as due to a 
't Hooft fermion, and it seems unphysical. In order to get more interesting and reasonable 



results, eliminating the poles at k 2 = in c[ nh ^ and Cj?^, we can take E\ — — 1 and 
E1/2 = —(1 + f P1/2), see eqs. (3.30) and (3.31). As opposed to eq. (4.7), the IR limit of 
the weighted sum —(Co — 4Ci/ 2 + 3Ci/ 2 ) depends now on four parameters, the p s and E , 
so that one can easily obtain a positive mass-squared for the sfermions. For defmiteness 



and for an easier comparison with eq. (4.7), consider taking p\ = 00 and p\li — ~ 2/3, so 
that 

2p + 2E 



- (Co - 4C 1/2 + 3d 



(4.9) 



>o 1 + 2p 

which can be positive if E < — p (assuming a positive po). The sfermion masses can 
then be even bigger than the Majorana gaugino mass if i/1/2 is somewhat smaller than 
V\Eo~\V- 

The punchline of the above analysis is that tuning appropriately the boundary con- 
ditions, one can realize holographically any scenario between pure gaugino mediation 



25 32 to minimal gauge mediation 33 35 as well as scenarios with suppressed gaugino 



masses 



36 42 which would fit into a split supersymmetry scenario 43,44 . Hence, hard 



wall models can actually cover all of GGM parameter space. In fact, it is not entirely 
satisfactory that a necessary ingredient for all this amounts to introduce two parameters, 
H\i% and E , which are directly proportional to gaugino and sfermions masses, respec- 
tively. This is reminiscent of minimal benchmark points. It would thus be desirable to 
try and obtain both Majorana gaugino masses and positive squared sfermions masses by 
enriching the dynamics in the bulk instead of introducing inhomogeneous terms in the IR 
boundary conditions. In the next section we will achieve this goal by turning on a linear 
profile for an R-charged scalar, as it was done in [I]. 

Let us finally mention that, as noticed in (45], a positive value for C\ — Co is a desirable 
feature, in that it helps raising the mass of the Higgs in gauge mediation scenarios. In 
our models, this is achieved by the same conditions which make the right hand side of 



(4.9) positive. 



5 Hard wall with R-symmetry breaking mode 

In this section we would like to construct a simple scenario in which the R-symmetry 
is broken (and gaugino masses generated) dynamically. We will follow the same logic 
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as [T], where it was observed that considering only the minimal action (2.1) for the vector 



supermultiplet it is impossible to break the R-symmetry by bulk dynamics, and get non- 
zero gaugino Majorana masses. As in the top-down model considered in [I], we will see 
that the dynamical breaking of R-symmetry implies automatically the absence of massless 
modes in Ci/ 2 . Notice that this physical consistency condition had instead to be imposed 
by hand, in the previous section. 

We introduce a new dynamical scalar field rj in the bulk with m 2 = —3, and treat it 
as a linear fluctuation around the hard wall metric. 

The action for rj at the linearized level is completely determined by its mass while 
the precise values of its couplings with the vector multiplet can be guessed by analogy 
with the N = 2 supergravity embedding considered in [I], based on the general results 
of (46|[47p 



n 2 r 

S ki „ = ^ J d 5 xVG(G^d, V d uV - 3 V 2 ) , (5.1) 

Sint = ^ J £x y — [( V + zd zV )(xX + XX) + zdtfM + &)] . (5.2) 

One might think that, in view of the possibility of constructing more general bottom- up 
models, it might be interesting to see what happens if we take arbitrary coefficients in the 
interactions term. On the other hand, asking for a gravity dual of a supersymmetric field 
theory (which then breaks supersymmetry spontaneously or by a soft deformation) puts 
severe constraints on the possible interactions. In fact, precisely the constraints dictated 
by supergravity. One can check that choices other than the interactions above do not give 
the right supersymmetric result in the deep UV. 

We demand the R-symmetry breaking mode i](z, x) to have a non-trivial profile in the 
vacuum which is independent on the boundary space-time directions in order to preserve 
Poincare invariance of the boundary theory. The most general solution to the resulting 
equations of motion for rj without k dependence is 

rj(z) = zr}o + z 3 fj 2 , (5.3) 

where r] Q and 772 are two arbitrary constants. These two constants can be fixed imposing, 
as usual, boundary conditions at z = and at the IR cut-off z = This strictly 

amounts to considering them as free parameters, which we will do in the following. 



7 With respect to |Tj we have rotated the spinors by a phase in order to have a real B 1 / 2 at the end. 
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Ci Ci 




Figure 1: The solid blue line is for r\ = 0, the dashed one on the left figure is for r\ = 0.1 z 3 and 
on the right figure for rj = 0.1 z. In these plots pi/ 2 = 3 and fj, = 1. Notice that turning on rj 
the massless pole disappears. 



The equations of motion for A are modified by the presence of the extra contribution 



(5.2 ) and become 



(zd z - §)% + zo'hi + zd z r])£ = , 

„ - . ( 5 - 4 ) 
zd z + §)£ + zcr^x + (»7 + ^77) x = • 



Comparing with eqs. ( 2.8| ), one can easily conclude that in the present case the leading 



boundary behavior of the fermionic field is not modified with respect to the minimal case 
(2.12). On the other hand, as discussed in [TJ, whenever r/o ^ 0, we have to modify the 
definition of the fermionic correlator denning Bi/ 2 according to 

(ja{k)jp{-k)} = - + 2rj e al3 , (5.5) 



while the expression for the non-chiral fermionic correlator (2.15) remains unchanged 



We now need to solve eqs. (5.4) by imposing (homogeneous) boundary conditions in 
the IR (that for generic pi/2 would give a massless pole when 77 = 0). Unfortunately, this 
cannot be done analytically, and we have to resort to numerics. Figures [T] and [2] contain 
our results. 

It is remarkable to see that when the R-symmetry is broken by a scalar profile, the 
pole in C1/2 disappears automatically. We note that the sfermion mass-squared is driven 
positive by the fact that C1/2 is still quite large near k — 0, at least as far as rj is a 
perturbation. If we stick to this model without playing with inhomogeneous boundary 
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Figure 2: On the left 77 = 0.1 z 3 , on the right rj = 0.1 z. In these plots pi/ 2 = 3 and fx = 1. 
When rj ^ a non-zero -B 1 / 2 is generated. 



conditions in the IR, it can be seen that we are able to explore a smaller region of parameter 
space. (Possibly, a larger portion of parameter space can be reached by playing with pi/2-) 
While the above analysis is done numerically, it would be nicer to have some analytical 
control on (at least) the low momenta behavior of the correlators, to see, for instance, 
how the pole in C\/2 disappears when the R-charged scalar is turned on. This analysis 
turns out to be possible if we also take the parameters 7] and 772 parametrically small, 
and we obtain 

^ i + §pi/ 2 v 
1 + 2P1/2 fc + ^ M V0 ,fj 2 

^ 1 + |Pi/2 8p 2 M^ 2 7) 



l + lpi/2 fc 2 + 4M, 



'/(l-'/2 



where 



M m52 =,„ + I -t|^^. (5.8) 

Notice that these formulae agree with the numerical plots in figures [T] and [2j The details 
of this computation are explained in appendix |Bj 

6 Conclusions 

In this paper we have considered holographic models of gauge mediation, and used hard 
wall backgrounds as a prototype to see whether and how strongly coupled hidden sectors 
can actually cover the GGM parameter space. 
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For a generic choice of boundary conditions at the IR wall, the resulting low energy 
spectrum is that of mediation scenarios with extra non-SM gauge sectors, where Z'-like 
gauge bosons acquire a mass due to symmetry breaking in the hidden sector, and mediate 
supersymmetry breaking effects to the SM. Tuning some parameters one can eliminate 
the composite massless modes emerging in the hidden sector recovering more standard 
gauge mediation scenarios, and cover all of GGM parameter space. 

Hard wall models can be seen as an effective way to mimic non-trivial bulk dynamics 
(for instance, several effects of a hard wall can be obtained with a non-trivial dilaton 
profile). In this sense, it would be interesting to see whether there exist more dynamical 
models where part of all GGM parameter space can be reproduced. The analysis of 
section 5 is an example of such a strategy. More ambitiously, it would be interesting to 
see whether relaxing the AAdS-ness of the background (which, strictly speaking, makes 
the dual hidden sector dynamics somewhat unrealistic), changes any of the results. 

On a more formal level, it might very well be that not all sets of physically distinct 
boundary conditions can find a counterpart in fully-fledged lOd string embeddings. This is 
a question which is very difficult to answer, but which certainly deserves further attention. 

We hope to return to these issues in the near future. 
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A Definition and properties of Bessel functions 

In this appendix we will briefly summarize the definitions and the basic properties of the 



Bessel functions following 48 . As we show in the paper, the equations for the fluctuations 
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of the fields in the vector multiplet can be recast in the Bessel differential equation 

(z 2 d 2 z + zd z - (z 2 k 2 + u 2 ))f(z, k) = (A.i) 

for v — 0, 1. Taking k 2 > 0, the two independent solutions of the above equation can be 
written as 

f(z, k) = Cl (k)I u (zk) + c 2 (k)K u (zk) . (A.2) 
The Bessel functions I v {x) and K v {x) can be defined as power series: 

r 2n+u 

n\2 v+2n (v + n)\ ' 



n=0 



VK ) 2 ^ n\2 2n - u 

n=0 

00 2n+u i i 

E n!2 w ( „ + n ). ('°g (f ) " + V - + » + !))• (A-*) 

n=0 v 7 

with 

^(n + l) = V--7, (A.5) 

where 7 ~ 0.58 is the Euler-Mascheroni constant. 

Some recurrence formulas involving I v (x) and K u (x) are the following: 

d x I (x) = h(x), 
d x K (x) = -K ± (x), 
xd x I\{x) + h(x) = xIo(x), 

xdxK^x) + K x {x) = -xK (x). (A.6) 

We now list the asymptotic behavior of the Bessel functions of interest. In the limit 
x> 1 an analysis based on the integral representation of the Bessel functions shows that 
at leading order 

'-(*> -SL • < A ' 7 > 



x— >oc 



x^-oc 



,2xJ 

independently on v. 



K u (x) ~ (^) 1/2 e-% (A.8) 
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In the limit a; < 1 the Bessel functions are well approximated by the first terms of 



their series expansion (A.3)-(A.4): 



Ux) ~ 1 + -x 2 + 



h(x) ~ -x + — x 3 + . . . , 



>o 2 16 



yx) l± - logx + log 2 - 7 + 



K (x 



11/ , 1 

- ~ + o x logx-log2- - + 7 ) +.... 

+o x 2 \ 2 



(A.9) 



(A.IO) 



B The IR limit of correlation functions 



In section 5 we have shown how the presence of a non-trivial profile for an R-charged scalar 
field, T], while providing a non-vanishing value for the R-breaking fermionic correlator 
B1/2, consistently removes the pole from the non-chiral fermionic correlator C\/2- The 
analysis was done by numerical methods. Here we show that one can actually study the 
IR behavior of holographic correlators analytically. 

We are interested in analyzing the behavior of the correlation functions for small k. 
More precisely, the relevant quantity is h/fi -C 1, where z = is the position of the IR 
wall, so that the limit can also be seen as moving the wall closer to the boundary. This 
suggests that if we just need to evaluate the behavior of the C s functions at low momenta, 



i.e. (3.15 )-(3.17), we can impose the IR boundary condition directly on the near-boundary 



expansion of the solutions, keeping only terms up to a mode high enough to match the 
order in k 2 at which we need the C s . Indeed, in previous sections we have seen that this 
limit is very easy to obtain when one has exact solutions, since it involves expanding 
Bessel functions near the origin, i.e. keeping only the near-boundary expansion. 

Let us illustrate this procedure with Co with homogenous IR boundary conditions. 



We just need to substitute the expansion (2.10) in the boundary conditions (3.8). We get 
1 



:(d log(A//i) + d ) + Po^(2d log(A//x) + d + 2d ) = 



that is 



Applying 



d = -d log(A//i) + 



Po 



2p 



t Sd 
'5d '' 



(B.l) 
(B.2) 
(B.3) 
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we obtain (3.15) right away and effortlessly. In order to reproduce eqs. (3.16) and (3.17) 



the only added difficulty is that we have to go one order higher in the expansion, if 
interested in both the 1/k 2 pole and the finite term. 

Notice that this procedure works because the equations of motion themselves are not 
modified with respect to the AdS ones. If we had O(fi) corrections to the metric (as in 
the example used in fil), it would be impossible to take small without introducing 
large corrections to the background metric and thus to the equations for the fluctuations. 

The case of an AdS hard wall with a scalar profile turned is a particular case. In order 
to prove that the pole in Cxii disappears when rj = r]QZ + fj2Z 3 is turned on, we should take 
the limit k — > in such a way to keep terms of the form (k 2 + t/q) -1 or (k 2 + Z^ -4 ^!) -1 - 
Therefore, the correct scaling is 



(B.4) 



and we should focus on the order e~ 2 in the small e expansion of C\/2- Keeping 7] small we 
also ensure that we can still use the AdS near boundary expansion for the fluctuations. 
The same kind of expansion can be done for the Bi/ 2 correlator, with the difference that 
it starts from the e -1 order. In both cases, the leading terms in the e expansion receive 
a non-trivial contribution both from r] ^ and from ?j 2 7^ and they are determined by 
keeping the near-boundary expansion 



y 3/2 



Z0 + J2&n + &n\og(zA))z 2n 



n=l 



(B.5) 



up to n = 1 and n = 2, respectively. The results for the order e -2 of C\/2 and the 
order e _1 in Bi/ 2 are reported in section 5. If one wants to go to the next order in e, 
which is order e° for C\/2 and order e for B\/ 2 , one should keep terms up to n = 3 in the 
near-boundary expansion. Let us stress that this e expansion is different from a simple 
expansion for small momenta. For instance, the finite k = term will receive contribution 
from arbitrary high orders in e, which in turn would require to keep arbitrary high terms 
in the near boundary expansion. Nevertheless, as long as rjo and fj 2 are kept small, the 



approximations (5.6)-(5.7) give a reliable information about the finite value at k — 0, as 



can be checked with the numerical results plotted in figures [T] and [2] 
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